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Abstract 

The purpose of the CODALEMA experiment, installed at the Nangay 
Radio Observatory (France), is to study the radio-detection of ultra-high 
energy cosmic rays in the energy range of 10^^ - 10^^ eV. Distributed 
over an area of 0.25 W, the original device uses in coincidence an array 
of particle detectors and an array of short antennas, with a centralized 
acquisition. A new analysis of the observable in energy for radio is pre- 
sented from this system, taking into account the geomagnetic effect. Since 
2011, a new array of radio-detectors, consisting of 60 stand-alone and self- 
triggered stations, is being deployed over an area of 1.5 km^ around the 
initial configuration. This new development leads to specific constraints to 
be discussed in term of recognition of cosmic rays and in term of analysis 
of wave- front. 
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1 Introduction 

The radio-detection of ultra-high energy cosmic rays (UHECR) may be consid- 
ered as one of the future tools, likely to supplement the traditional measures of 
extensive air showers (EAS) [Ij. This method exploits, through antennas, the 
radio signal that accompanies the development of the EAS. Many results have 
highlighted its feasibility and its potential in reconstructing several parameters, 
as the arrival direction, the impact location, or the radio lateral distribution 
function of the electric field (RLDF) [2l|3l[4]. However, except a recent study [5j, 
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the estimation of the energy through radio has been only briefly discussed [6l[7], 
although this step is unavoidable if one really plans to use this technique. Simi- 
larly, the radio wavefront characteristic is remained poorly determined |8j[9l[T0], 
although its knowledge could be consider as one of the flrst steps in retrieving 
information about the EAS itself. These issues motivate the development and 
the continuation of demonstrative experiences, such as CODALEMA [llj. 
This paper summarizes some new aspects, especially regarding the energy esti- 
mator [12j and the point source of the radio emission [13j. 

2 Experimental setups 

The study of the energy response in radio took advantage of the original setup 
which uses a centralized trigger (see [14J). It consists of two arrays of detectors 
(flg. [T]): an antenna array of 24 short active dipoles distributed on cross of 
600 X 500 m; a particle detector array of 17 plastic scintillator plates spread over 
a square grid of 340 x 340 m. The antennas operate in the frequency bandwidth 
1 — 200 MHz. The scintillating detectors are read by 2 photomultipliers (one 
at high gain, the other at low gain) so as to cover an overall dynamics of 0.3 to 
3000 vertical equivalent muons (VEM). 




Figure 1: Aerial view of the CODALEMA setup installed at Nangay. On the 
original system, the plastic scintillators are represented by blue circles while 
yellow and orange circles represent the dipole antennas, placed respectively in 
the EW and NS directions. The white squares show the locations of the au- 
tonomous stations of the new array in deployment (33 are currently in opera- 
tion). 

The trigger system is provided by the coincident detection of the flve cen- 
tral particle detectors, leading to a detection threshold around 10^^"^ eV. For 
each trigger, all the sensors signals are digitized at a sampling rate of 1 GS/s 
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(snapshot of 2.56 /as). 

During the off-hne analysis, the radio transient is searched in the radio wave- 
form through the 24 — 82 MHz bandwidth, from which the electric field peak 
value and the arrival time are extracted. If more than 3 radio transients are 
detected in coincidence, the arrival direction of the electric field (zenith angle ^, 
azimuth angle (j)) is reconstructed by considering a planar wavefront. Assuming 
a RLDF for the electric field of the form sq ex.p{—d/do) [15], four additional 
observables are derived: the electric field magnitude at the radio shower center 
£0, the slope parameter of the electric field do, and the impact location of the 
radio shower center on the ground (xq, ^o)- The standard deviation of the noise 
on each antenna is used as the error on the electric field measured. From Monte- 
Carlo simulations, an mean error of 22 % is estimated for eq. The timing error 
is evaluated to 10 ns, leading to an angular accuracy of reconstructed arrival 
directions in the range of the degree. 

Meanwhile, the arrival direction of the particle front (assumed planar) is com- 
puted from the relative arrival time of each particle detector and an analytical 
NKG lateral distribution [T6] is adjusted on the measured particle densities in 
the shower frame. If the shower core position falls inside the particle detector 
array, the event is referred to as "internal" and a adjusted NKG distribution 
is used to deduce the primary energy using the constant intensity cut (GIG) 
method [T7J . This processing leads to a "particle" energy resolution of the order 
of 30 % at 10^^ eV. 

A coincidence below 100 ns in time and 20° in arrival direction with the event 
"particles" labels the event "radio" as an EAS. Imposing zenith angles less than 
60°, and energies "particle" above 10^^"^ eV, 376 coincidences were available for 
studies of the energy "radio", discussed in section 3. 

The study of the wavefront features has been based on the new autonomous 
setup made of self-triggered antenna stations (see |18j). The later are based 
both on the butterfy antenna and on a embedded electronic. Its aim is to 
complete the initial configuration, providing finer measurements of the electric 
field (with two polarization states per station: EW and NS), better SNR, and 
an increased range in energy, using a larger area. 33 stations are currently in 
operation around the original setup. Ultimately, this apparatus will cover an 
area of 1.5 /cm^, with the additional aim to bring this technique reliable and 
mature for large infrastructure. Up to now, efforts are mostly focused on the 
radio- sources, observed using the self-triggering technique. The interpretations 
of these observations are central to the fourth section. The event are built from 
the time coincidences between the stations using the GPS time. 

3 Energy estimation 

This study has focused on the possibility of deriving from radio, a simple ob- 
servable which could provide a preliminary satisfactory estimate of the energy of 
the shower, by going further in some previous works [7, 14j. Because the electric 
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Figure 2: Electric field Sq versus the primary energy Ep (for readability, the 
error bars of 30 % for Ep and 22 % for eq are not represented). The result of 
the adjustment of the data with the power law Ep^^ has motivated the use of 
a linear adjustment depicted in red. 

field radiated should be linked to the energy of the primary particle through the 
number of secondary charges produced within the shower, we chose to exploit 
the maximum electric field Sq (derived from the RLDF) as radio-observable of 
the energy. For this, the correlation with the shower energy estimator Ep (de- 
duced from the CIC analysis) is exploited (fig. [2| by fitting the scatter plot 
distribution {Ep^ sq) with the linear relation sq = aEp + f3. 

Inversion of the latter gives the energy estimator of the radio Eq = eo/a — 
P/a. Discrepancies between both observables Eq and Ep is quantified through 
the analysis of the distribution of their relative differences {Ep — Eq)/Ep char- 
acterized by its standard deviation (SD). A typical forms of the distribution 
{Ep — Eo)/Ep is presented in fig. [sj 

To accord with our earlier works [19j, which showed that the electric field 
magnitude depends on the Lorenz force via the vector cross product v A B, we 
have corrected the estimator by the projection of the Lorentz force along the 
EW polarization. So, two detected showers of the same actual primary energy 
E (and of the same nature) will produce two different radio signals depending on 
their angle with the geomagnetic axis, by following the relation: Sq = E/\{v A 
'B)ew\' In doing so, the scatter plot (£^p, €o/\{'v A 'B)ew\) (fig- |4| shows that 
events located close to the earth's magnetic field direction (|(v A B)^;^^! < 0.1) 
gain unrealistic energies, while their CIC energies show nothing unusual. The 
original distribution of (|(v A B)^;^^^! versus Ep is presented figjs] 

To account for this, we hypothesized that the observed electric field Sq can 
also be proportional to the actual shower energy E (to produce observable sig- 
nals, especially when the Lorentz contribution becomes weak). This leads to 
modify eo toward Sq 6:o/(|(v A B) ew\ + c) where c is a scalar quantity. With 
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Figure 3: Distribution {Ep — Eo)/Ep obtained after calibration and for the full 
sample of events (blue crosses). The standard deviation is 31 % with c = 0.95 
(see below). The red line shows a Gaussian fit made in the interval [—1, 1]. 




Figure 4: Scatter plot of the corrected electric field £:o/|(v A B)^;^/^! versus the 
primary energy Ep . The linear adjustment is shown in red. The green arrows 
show some of the typical over-corrected events. The systematic analysis of their 
features indicates that all these events have arrival directions near the earth's 
magnetic field direction (|(v A B)£;v[/| < 0.1). (These points are removed for the 
linear fit of the scatter plot). 
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Figure 5: Distribution of |(v AB)^;^!/! versus the primary energy Ep for the raw 
sample of events. 



the quantity |(vAB)£;^| ranging from to 1, c = corresponds to a purely geo- 
magnetic contribution, while a high c value (e.g. c > 6) should represent the case 
for which the geomagnetic emission becomes negligible compared to the second 
effect. A high c value must induce a behavior similar to that observed without 
correction. The figure [g] presents the standard deviation of {Ep — Eq)/Ep as a 
function of c and for several bins in |(v A B)£;vt^|. 

The observed variations in function of |(v A B) ew \ + c suggest strongly that 
the radio signal is made of the mixture of several emission effects. Because the 
observable eo seemed sensitive enough to different contributions, its association 
to an existing electric field was then studied, for instance to an electric field 
component whose orientation would be linked to the radio shower axis. To 
address this possibility, we changed the term c by: £q ^ £^|(v A B)£;vt/| + 
E c\ sinO sin^|. The effect of this correction is depicted fig. [7| which presents 
the standard deviation of {Ep — Eq)/ Ep as a function of c and for the same 
windows as in |(v AB)£;v[7|. A worse standard deviation is obtained, which leads 
to exclude the scenario of a longitudinal contribution, relative to the shower axis. 

Some recent developments [H |20l [21] could however give room to the simple 
correction with c. For instance the charge excess mechanism could be in qual- 
itative agreement with the presented effects. Indeed, this emission depends on 
the charges produced in the shower (and so from the energy) and its weight in 
the total emission process should increase when the geomagnetic contribution 
decreases. The effect of variation of the optical index is also expected to make a 
contribution rather independent of the geomagnetic angle, although our data do 
not suggest this contribution in the RLDFs. On another view, our result could 
also suggest a coherence effect modulated by the Lorentz force. The coherence 
mechanism should increase when the secondary particles are less dispersed by 
the Lorentz force (i.e. at small |(v A B)£;v^|) during their propagation. 
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Figure 6: Evolution of the Standard Deviation of the distribution {Ep — Eo)/Ep 
with the correction factor 1/(| (v A B) ew I + c) as a function of c and for several 
windows of |(v A B)^;^^/^!. The blue line represents the whole sample of events. 
The stars located at c = 10 indicate the values of the SD without correction. 
The SD errors are given by V2 crVVATTI (TV ^ 70 for the subsets). 
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Figure 7: Evolution of the Standard Deviation of the distribution {Eq — Ep)/ Ep 
with the correction factor l./(|(vAB)£;vt^|+c| sin sin (/) | ) as a function of c and 
for several values of |(v A B)£;vf|- The blue line represents the whole sample of 
events. The stars located at c = 10 give the value of the SD without correction. 
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4 Source reconstruction 

Due to its potential link with the quantity Xmax^ the reconstruction of an appar- 
ent point source in the emission could be one of the significant key-point in the 
development of radio-detection [8j. Because the arrival timing on each antenna 
is determined at the maximum amplitude of the radio signal, itself dominated 
by the maximum of the shower [22j, this observable should therefore be linked 
to this part of the longitudinal development of the shower. The assumption of 
a spherical wavefront, then gives, in a natural way, a efficient tool to find the 
position of this instant of emission. Moreover, the self-trigger technique in radio 
can be subject to delicate limitations in regard to UHECR recognition, due to 
noises source interferences (human activities, stormy conditions). In this con- 
text, the identification of the positions of the sources of emission is crucial. For 
sources outside of the setup, this identification can be only achieved by using 
curved wavefront. 

For a spherical wavefront, the commonly used technique to extract this 
curvature radius relies on the calculation of the function, using the ar- 
rival times and the locations of the antennas, by minimizing the quantity: 
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X^i^l^ts) = \\rs - TiW - {ts - U) with {ri.U) the location and re- 
ception time of the i^^ antenna; {jl^tg) the searched location and unknown 
instant of emission of the emitting source; and c the wave velocity consid- 
ered constant. Figure [S] shows a typical reconstruction of sources obtained 
with CODALEMA, by invoking such a minimization of a spherical wave with a 
Levenberg-Marquardt (LVM) algorithm. 

The source positions are poorly reconstructed, although they are fixed and 
that the number of events is high. Additional simulations showed that what- 
ever the samples (versus any distances and arrival directions of source, time 
resolutions, initial conditions, and detector configurations), different algorithms 
(LVM, Simplex, Line-search) lead to similar conclusions (see fig. [9|. It suggests 
this functional seem ill-adapted to the common minimization algorithms. 

Applied to a cosmic event, which is a single realization of the observables 
detected (arrival times and peak amplitudes of antennas), such drawbacks could 
still worsen, even by extracting mean values by using statistical approaches. To 
understand the observations, we have undertaken to study the main features of 
this function [23j. 

Using the criterion of Sylvester [25] and the analysis of the principal minors 
of the Hessian matrix, we found that the x^ function is not convex on small 
domains, and thus is likely to exhibit several local minima, according to source 
position relative to the antenna location (fig. [1Q| ). These minima induce prob- 
lems to converge to the correct solution. 



The development of the first optimality condition (Jacobian = 0) allowed to 
study the critical points. It indicates that the iso-barycenter of the antenna ar- 
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Figure 8: Typical result of reconstruction of two entropic emitters at ground, 
observed with the stand-alone stations of CODALEMA (red circles), through a 
standard minimization algorithm. Despite the spreading of the positions recon- 
structed, both actual transmitters are stationary point sources. Such patterns 
are also commonly observed in others radio experiments [24| . 




Figure 9: Results of the reconstruction, using a testing array of 5 antennas 
(400 X 170 m), of a source with a radius of curvature of 10 km and the Simplex 
algorithm. The blue line shows the distribution of the initial values for the blind 
search. 
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Figure 10: Plotting of the objective-function versus |r^| in the case the same 
testing array as fig. [9] for a source outside the convex huh. This configuration 
leads to multiple local minima, located on a line joining the antenna barycenter 
to the true source. 

ray (of the lit antennas for a given event) plays an important role in explaining 
the observed numerical degeneration. There are strong indications, in agree- 
ment with the experimental results and our simulations, that the critical points 
are distributed on a line connecting this barycenter and the actual source loca- 
tion. In this case the problem is ill-posed. 

This approach has highlighted the role of the convex hull of the array of an- 
tennas, which is defined, event by event, by the envelope of the tagged antenna 
(fig. [TT]). We found that if the source is in the convex hull of the detector, the 
solution is unique. In contrast, the location of the source outside the convex 
hull of the detector, causes degeneration of solutions (multiple local minima) 
regarding to the constrained optimization problem. 

We have used this observations to construct an alternative method of locat- 
ing the source. To avoid the trap of the local minima with common algorithms, 
we chose to compute directly the values of the objective function on a grid, using 
a subset of phase space in the vicinity of the solution a priori, and assuming that 
the minimum of the objective function corresponds to the best estimate of the 
position of the source of emission. The direction-priori is given by the planar fit, 
while the quantity is left free in the range 0.1 — 20 km (the upper bound being 
determined by the value of the curvature exploitable, given the time resolutions 
currently available). A typical result obtained with our method is presented in 
fig. [12] The first results indicate better accuracies. Next step of this study 
will consist to exploit all the available data (events other than "internal"), and 
quantify the bias introduced by this method. 
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Figure 11: Sketch of the topology of a convex huh of a 2D array of antennas for 
a particular event. The dark markers are the lit antennas for this event. The 
red dot line shows the corresponding convex envelope. The red cross depicts 
the centroid of the lit antennas, the black cross represents the barycenter of the 
antenna array. 
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Figure 12: Histogram of the reconstruction of a source located at 10 km from 
the same the detector array as in fig. 8, and using the grid method associated 
to the search of the global minimum for each event. 
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5 Conclusion 

The possibility to extract in radio signal, a simple observable providing an esti- 
mate of the energy of the shower has been presented. The maximum amplitude 
of the radio electric field on the ground (eo) seems relevant (and surprisingly very 
usable when compared to particles detectors). We show that this observable is 
sensible to different mechanisms causing the radio signal. We demonstrate that 
the only correction by the geomagnetic effect would lead to an over-evaluation 
of the energy of the EAS. We deduce that a further contribution is also at the 
origin of the formation of the radio signal. In a first interpretation, effects like 
coherence, charge excess, or Cerenkov emission could play a significant role in 
the observed amplitude modulation. More data should improve obviously the 
proposed approach, but given the results, the radio signal may be sufficiently 
sensitive to the longitudinal evolution of the shower, to contain valuable infor- 
mation on the nature of the primary. 

In second part, we addressed the problem of the apparent point of emission in 
radio, by reconstructing spherical wavefronts. Indeed, the experimental results 
indicated that the usual methods of minimization of spherical wavefronts could 
induce a mis-localization of the emission sources. The specific characteristics 
of the function can lead to the formulation of an ill-posed problem, in the 
sense of minimization theories under constraints. The conventional algorithms 
can exhibit a degeneration of the solutions, due to the presence of several local 
minima. To avoid possible traps in these minima, we tried a simple method. It 
is based on a grid calculation of the function ^ind the search for an absolute 
minimum in a phase space fixed a priori. Further developments and tests will 
undoubtedly still needed to use the method routinely. 
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Pays de la Loire", under grant 2007/10135. 
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